The Tolman-Bondi-Vaidya Spacetime: matching timelike dust to null dust 
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The Tolman-Bondi and Vaidya solutions are two solutions to Einstein equations which describe 
dust particles and null fluid, respectively. We show that it is possible to match the two solutions in 
one single spacetime, the Tolman-Bondi-Vaidya spacetime. The new spacetime is divided by a null 
surface with Tolman-Bondi dust on one side and Vaidya fluid on the other side. The differentiability 
of the spacetime is discussed. By constructing a specific solution, we show that the metric across 
the null surface can be at least C 1 and the stress-energy tensor is continuous. 
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INTRODUCTION 



The Tolman-Bondi solution describes the inhomoge- 
neous gravitational collapse of spherical dust. In comov- 
ing coordinates, the metric is given by 0,0 



ds 2 



-dt 2 



1 + E(R) 



dR 2 +r 2 dtl 2 , 



(1) 



where d£l 2 is the metric of the unit two-sphere, E(R) is an 
arbitrary function of the comoving radial coordinate R, 
and the prime denotes the partial derivative with respect 
to R. The function r(t,R) is a solution of 



2M{R) 



E(R), 



(2) 



where the overdot denotes the partial derivative with re- 
spect to t and M(R) is the effective gravitational mass 
within R. In contrast, the incoming Vaidya metric de- 
scribing collapsing null dust is given by 0,3 



ds 2 



2M(v) 



dv 2 + 2dvdr + r 2 dn 2 



(3) 



where v is an advanced time and the v — constant line 
gives the trajectory of the null fluid. Both the Tolman 
and the Vaidya spacetimes provide important examples 
for the study of naked singularities formed in gravita- 
tional collapse. The intrinsic relation between the two so- 
lutions had not been investigated for many years. Lemos 
first showed in [j| , that the Vaidya metric © can be ob- 
tained from the Tolman-Bondi metric Q by taking the 
limit E — > oo. Therefore, the two metrics belong to one 
family and their naked singularities are of the same na- 
ture. Hellaby [6] pointed out that Lemos' derivation re- 
quires that E be a constant through the spacetime. As a 
generalization, Hellaby showed that this assumption can 
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be dropped and the same Vaidya limit can be achieved 
from the Tolman-Bondi solutions by using a different co- 
ordinate transformation. 

Despite the connection discovered above, the two solu- 
tions have been treated in two separate spacetimes which 
are connected only mathematically by taking the limit on 
E. Is it possible to put the two solutions in one single 
spacetime such that the Tolman-Bondi dust can approach 
the Vaidya fluid as a physical limit? Hellaby's work 
has provided some necessary mathematical framework to 
solve this problem. Hellaby makes no assumptions on 
the function E{R) in the Tolman-Bondi model. This is 
an important point for our generalization. In the 
Vaidya limit is obtained by taking the limit E(R) — > oo 
everywhere. 

Our approach to matching the two solutions in one 
spacetime is simple. Instead of requiring E{R) — > co 
everywhere, we require E(R) to diverge only at one hy- 
persurface R — Rq. For R < Rq, E(R) is finite and it cor- 
responds to a Tolman-Bondi solution. For R> Rq, the 
spacetime transfers into the Vaidya null fluid. In this new 
spacetime, dust particles which follow timelike geodesies 
and photons which follow null geodesies coexist and a 
phase transition occurs across the R = Rq hypersurfacc. 
By imposing some general conditions, the spacetime is 
continuous (C°) by construction. An important issue to 
be discussed is the smoothness around the limit surface 
R = Rq, i.e., how well one can match the two solutions 
together. We show from a specific example, that the met- 
ric across the limit surface can be at least C . We also 
show that the stress-energy tensor and the Kretschmann 
scalar are both continuous in the spacetime. 



II. CONSTRUCTION OF THE SPACETIME 

In this section, we shall construct a spherically sym- 
metric spacetime which combines the Tolman-Bondi dust 
and the Vaidya null dust. To be definitive, the Tolman- 
Bondi dust will be configured inside the Vaidya dust. 
However, all of our arguments and results apply to 
the reversed configuration. We start with the Tolman- 
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Bondi metric (JTJ. It is evident from this metric that 
each comoving observer follows a timelike worldline with 
R — constant. To approach the Vaidya limit, we shall let 
E{R) diverge at some fixed R — Ro > such that the 
spacetime consists of the Tolman-Bondi dust for R < Ro 
and the Vaidya dust for R > Rq. If this can be done, 
the R = Rq surface will become a limit null surface for 
Tolman-Bondi observers, which indicates that the origi- 
nal Tolman-Bondi coordinates in Eq. (JTJ will break down 
at R = Rq. This is reminiscent of the Schwarzschild solu- 
tion where the original Schwarzschild coordinates break 
down when one approaches the event horizon, and the in- 
troduction of the Kruskal extension is necessary to elim- 
inate the apparent coordinate singularity. Now we face a 
similar situation and we shall introduce a new set of coor- 
dinates to extend the spacetime from the Tolman-Bondi 
interior to the Vaidya exterior, to yield a Tolman-Bondi- 
Vaidya spacetime. Following 0, we rewrite Eq. 0J, 
eliminating dt through 

; , , —r'dR + dr , 
dr = fdt + r'dR -> dt = : (4) 



where a = a(R) is another arbitrary function. We shall 
show in section IlII Bl that t can be written as an explicit 
function of (R, r). Thus, each component in the metric 
of Eq. J5J depends on (R, r, 8, <f>). Hellaby required that 
M and r be finite while E goes to infinity everywhere. 
Here we adopt a similar requirement with only a small 
modification, albeit important: we require M and r re- 
main finite everywhere and E(R) is divergent for R > Rq. 
Then we have the following limiting forms at R = Rq 



'E(a - t) 

-Ve 

rE' , r= 
2E +a ^- 



(7) 



By simple substitution, we find the following limits for 
the components of Eq. JSJ (see |||) 



1 



and then using Eq. @, to obtain 



ds 1 



2M\ r' 2 
~r~) [1 + E{R)\ r 



-dR 2 + 2—dRdr 



-dr 2 + r 2 dn 2 



(5) 



Now the metric is expressed in terms of (R, r, 8, 4>). Since 
r is the areal radius, it should be kept as a good coor- 
dinate in the extended spacetime. We note, in passing, 
that Hellaby further introduced a new coordinate v de- 
fined by v = J Q R -j^dR — , where a is some function 
of R, such that the metric is finally expressed in coordi- 
nates (v, r, <fi). However, the components of the metric 
involve functions like E(R) and M(R) that are implicit 
functions of (v,r). Consequently, it will be difficult to 
study the differentiability of the metric with Hellaby's 
coordinate system. Indeed we find that introducing the 
new coordinate v is not necessary here. Under certain 
conditions, coordinates (R, r, 8, <ft) are good enough for 
our purposes to describe the spacetime. To proceed, we 
first give some limiting forms of the components of the 
metric (jSJ) which has been derived by Hellaby in jfj]. 
The exact solution for r' in a collapsing model is Q 



'AT 
~M 



~E 



M 3E' . 

)(a-t) 

M 2E ' y ' 



(6) 



iL f^L rE ' \ 

t 2 WE + 2E 2 ) ' 



(8) 



E 



a' rE'Y 

(1 + E)r 2 \1 + EJWe + 2E 2 ) 



If we further impose restrictions such that 



a' rE' 
71 + 



1 at R — Rq , 



(9) 



and identical to 1 for R > Rq, one sees immediately, 
by comparing the components of the original Tolman- 
Bondi metric JSJ and those of the Vaidya metric © and 
identifying v with R, that the metric of Eq. (JSJ becomes 
exactly the Vaidya metric for R > Ro. In summary, 
our spacetime is constructed in the following way: The 
coordinates are (R,r,8,<fr). For R < R , the spacetime 
is described by any Tolman-Bondi dust in the form of 
(J5J which satisfies the limiting condition J^J and E — > oo 
at R = Rq. For R > Rq, the spacetime is described by 
the Vaidya null dust © with v relabeled by R. In this 
case, we can choose a smooth mass function M(R) for 
the entire spacetime. With this choice, the spacetime 
metric is continuous (or C°) by construction. For future 
reference, we write the complete metric as follows 



f x « _£_ dR 2 + 2 r_ dRdr _ ^2 + r 2 d tf R < Rq 

' ] L - 2 -MM) dR 2 + 2dRdr + r 2 dCl 2 , R > R . [ > 



r 



We see that the comoving coordinate R of the Tolman- 



Bondi dust automatically becomes the comoving coordi- 
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nate of the Vaidya null dust. This shows the unification 
of the two solutions expressed in this coordinate system. 
It is also easy to check that the RR component of the 
inverse metric, g RR , approaches zero as R — ► Ro, which 
confirms that the hypersurface R — Rq is a null hyper- 
surface. In the next section, we shall show through a 
concrete example that the metric at R = Rq can be at 
least C 1 . 



rior (R < Rq) and be identical zero in the Vaidya exte- 
rior. Hence, we choose h(R) as 



h{R) 



(Ro 




R) 3 R<R 

R > Rq. 



(12) 



To satisfy the limit we define a(R) in the Tolman- 
Bondi interior as 



III. AN EXAMPLE 

To demonstrate how well the two metrics can be 
matched, we provide an explicit example based on the 
prescription given in the previous section. 

A. Choosing the arbitrary functions 



a(R) 



\/Rq — R 



(13) 



In the following two subsections, we shall investigate the 
differentiability of the metric associated with the model 
spacetime and the continuity of the Einstein tensor across 
the null surface E. 



Since our interest is matching the Tolman-Bondi and 
Vaidya solutions, we shall restrict our attention to the 
neighborhood of the limit null hypersurface, which is de- 
noted by E. Since at R — Ro the function E(R) goes 
to infinity, it is convenient to discuss the smoothness of 
E(R) at R = Rq through an auxiliary function h(R) 
given by 



h(R) = 



1 



E(RY 



(11) 



According to the discussion in the previous section, h(R) 
should approach zero at E from the Tolman-Bondi inte- 



B. Differentiability of the metric 

Our task in this subsection is to show that the metric 
is C 1 on the null surface E. We shall show that the 
components of the metric is differentiable on E. 

We begin with the RR component in Eq. (|10|l . de- 
noting it on the Tolman-Bondi side, R < Rq, by g RR , 
and on the Vaidya side, R > Rq, by g RR - Combining 
Eqs. {IU, tEl an d on e may write g RR as 

a function of (R, r) 



2M(R.) 



9rr 



( 1 4- 1 ,.. "> f 1 i 2M(R) \ 



-t(R,r) 



(R -R) 3 / 2 \VRo~R 



3 M'(R) \ 
R-R Q + M(R) J 



M 



2R - 2Rq ' M(R) J 



(Ro - R) 3 



-I 2 



2M(R) 



(14) 



r 



where we have expressed t formally as a function of (R, r) 
To find the explicit form of i(i?,r), we consider the hy- 
perbolic solution for the collapsing model 6] 

^-(cosh?7 — 1) , 



(sinh T] — rj) = 
Solving Eq. I|15l) for rj gives 



E 
E 3 ' 2 (a 



f) 



cosh 



1 



M 
Er 



(15) 
(16) 

(17) 



where cosh is the inverse function of cosh and is re- 
quired to be positive. By symmetry, there is also a neg- 



ative solution for 77. But one can easily check that, as 
E —> 00, the negative rj leads tor^ \ff(t — a), which is 
not a collapsing solution. Therefore 77 is large and posi- 
tive in the region of interest, which implies that sinhr; is 
positive too. Thus, from Eqs. 115|l and (|16fl . we find 



JEr(2M + Er) - Mcosh -1 (l + Er/M) . 
t = a ^3f 2 . (18) 

It is obvious that t diverges at R = Rq for the functions 
chosen in section IlII Al This confirms the breakdown of 
the original coordinates, as we have mentioned at the 
beginning of section [H] Now we have g RR as an explicit 
function of (R,r). Its derivatives can then be calculated 
straightforwardly. To obtain the limits of these deriva- 
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tives at R — Rq, we treat Rq — R as a small quantity 
and compute the power series expansion about R = Rq. 
This treatment guarantees the accuracy of the final re- 
sults because no approximation has to be made in the 
intermediate steps. With the help of Mathematica, the 
results can be obtained almost instantly. We find that at 

R = Rq, 



9 9rr 


d 9 RR 


2M(R ) 


dr 


dr 




d 9 RR 


9 9rr 


2M'(Ro) 


dR 


dR 


r 



(19) 



Therefore, the RR component of the metric is differen- 
tiable (or C 1 ) at E. 

With similar notations, we repeat the procedure above 
for g^ r and g~ r and find that their derivatives with re- 
spect to R and r vanish as R — > Rq, agreeing with their 
Vaidya counterparts (note that g^ r and g+ r are both con- 
stants). Therefore, we have shown that the spacetime 
metric is C 1 across the null hypersurface E. 

However, it is worth mention, as a remark, that there 
are examples where the metric is not always C 1 on E 
when the conditions in section [H] are satisfied. Here is 
a counterexample. We choose, in the R < Ro region, 
h(R) = (Rq - R) 2 and a(R) = - ln(Ro -R). It is easy to 
check that the limiting condition (J5J is satisfied for this 
choice. By repeating the similar calculation we carried 
out for the first example, we find that dg^ r /dR, — ► — r 
when R — > R , which is not equal to the g^ r /8R — > 
limit except for the r = singularity. Thus, this example 
shows that the metric is not differentiable (or C 1 ) on 
E. There is a qualitative explanation for why the first 
example works while the second one does not. Since h(R) 
vanishes identically for R > Rq, we see h(R) defined by 
Eq. (O is C 2 at R = Rq while h(R) defined above is 
just C 1 . Roughly speaking, a better behaved function 
corresponds to a better behaved metric. 



Continuity of the stress-energy tensor and the 
Kretschmann scalar 



We have shown that the spacetime metric defined at 
the beginning of section [H] is at least C 1 . In this sub- 
section, we shall show that the spacetime performs even 
better than C 1 . One important quantity associated with 
the second-order derivatives of the metric is the stress- 
energy tensor T a b, which can be derived from Einstein's 
equation 



Gab — Rab ~ ^Rg a b — &nT ab 



(20) 



We now calculate G a b associated with the metric defined 
by Eq. H1(J|) and functions chosen in section UlI Al From 
Einstein's equation, it is straightforward to check that the 
only nonzero component of G a b for the Vaidya metric is 



We need to show that G a b associated with the Tolman- 
Bondi metric agrees with this limit at R = Rq. The cal- 
culation follows the same strategy as in subsection IIIIBI 
The limits are taken only in the final step by means 
of power expansion. We finally find that the Tolman- 
Bondi stress-energy tensor approaches the Vaidya limit 
at R = Rq, i.e., all components of G a b vanish at R = Rq 
except G R R^2M'(Ro)/r 2 . 

Another quantity related to the spacetime curvature is 
the Kretschmann scalar K = R abcd R a b c d, which is given 
in |6( for the Tolman-Bondi metric 



T 48M 2 32MM' 12M' 2 , s 

K T = — — 1 -r-^ (22) 



and for the Vaidya metric 



48A/ 2 



(23) 



G 



RR 



2M'(R) 



(21) 



Since r' — > oo in our model spacetime, we see that the 
Kretschmann scalar is continuous across E. 



IV. CONCLUDING REMARKS 

We have constructed, by generalizing Lemos and 
Hellaby's methods 0, , a spacetime where the Tolman- 
Bondi dust and the Vaidya null fluid coexist in harmony. 
The spacetime is described by a uniform coordinate sys- 
tem which is comoving with both the Tolman-Bondi and 
Vaidya observers. We have shown, from an explicit ex- 
ample, that the spacetime metric is at least C 1 across 
the null surface E. Moreover, we have shown that the 
stress-energy tensor and the Kretschmann are continu- 
ous on E. This property fails for some other spacetimes. 
For instance, when matching the Friedmann interior to 
the Schwarzschild exterior p|, the energy density across 
the matching surface is obviously discontinuous since the 
Friedmann interior has a uniform distribution of density 
and the Schwarzschild exterior has a vanishing density. 
Although the Tolman-Bondi and Vaidya models are two 
well-known solutions to Einstein's equations, the match 
of the two models is a new issue. We have studied the 
possibility of matching them in one spacetime. It re- 
mains unknown to us how this match can be realized in 
a physical process, or what physical mechanism can make 
it happen. 
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